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We present the calculation of the decay width of the Z-boson into three jets including complete quark mass 
effects to second order in the strong coupling constant. The study is done for different jet clustering algorithms 
such as EM, JADE, E and DURHAM. Because three-jet observables are very sensitive to the quark mass we 
consider the possibility of extracting the bottom quark mass from LEP data. 



1. INTRODUCTION 

A precise theoretical framework is needed for 
the study of the quark mass effects in physical 
observables because quarks are not free particles. 
In fact, the quark masses should be seen more 
like coupling constants than like physical param- 
eters. The perturbative pole mass and the run- 
ning mass are the two most commonly used quark 
mass definitions. The perturbative pole mass, 
M(p 2 = Af 2 ), is defined as the pole of the renor- 
malized quark propagator in a strictly perturba- 
tive sense. It is gauge invariant and scheme in- 
dependent. However, it appears to be ambiguous 
due to non-perturbative renormalons. The run- 
ning mass, fh(fi), the renormalized mass in the 
MS scheme, does not suffer from this ambigu- 
ity. Both quark mass definitions can be related 
perturbatively through 



M = fh(n) \ 1 + 



<*s(m) 



log 



m 2 (/i) 



(1) 



Heavy quark masses, like the bottom quark 
mass, can be extracted using QCD Sum Rules 
or lattice calculations from the quarkonia spec- 
trum, see [0 and references therein. The bot- 
tom quark perturbative pole mass appears to 
be around M b = 4.6 — A.7(GeV) whereas the 
running mass at the running mass scale reads 
fh b (fh b ) = (4.33 ± 0.06)GeU. Performing the 
running until the Z-boson mass scale we find 
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fh b (M z ) = (3.00 ± 0.12)GeU. Since for the bot- 
tom quark the difference between the perturba- 
tive pole mass and the running mass at the Mz 
scale is quite significant it is crucial to specify 
in any theoretical perturbative prediction at Mz 
which mass should we use. 

The relative uncertainty in the strong coupling 
constant decreases in the running from low to 
high energies as the ratio of the strong coupling 
constants at both scales. On the contrary, if we 
perform the quark mass running with the extreme 
mass and strong coupling constant values and 
take the maximum difference as the propagated 
error, induced by the strong coupling constant 
error, the quark mass uncertainty increases fol- 
lowing 



e r (m(M z )) ~ £ r (m(/u)) 



27o ( «s(m) 



(2) 



Po \as(M z ) 



- 1 e r (a s (M z )), 



where 70 = 2 and fto = 11 — 2/3Np, see figure [I]. 
We use the world average || value as = 0.118 ± 
0.006 for the strong coupling constant. 

It is interesting to stress, looking at figure [l], 
that even a big uncertainty in a possible evalua- 
tion of the bottom quark mass at the Mz scale 
can be competitive with low energy QCD Sum 
Rules and lattice calculations with smaller er- 
rors Furthermore, non-perturbative contribu- 

2 A recent lattice evaluation jlo| has enlarged the initial 
estimated error on the bottom quark mass, 7715(7715) = 
(4.15 ± 0.20)GeV, due to unknown higher orders in the 
perturbative matching of the HQET to the full theory. 
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tions are expected to be negligible at the if-boson 
mass scale. 

The running mass holds another remarkable 
feature. Total cross sections can exhibit poten- 
tially dangerous terms of the type M 2 log M 2 Is 



k 

that however can be absorbed |3j using Eq. ([y) 
and expressing the total result in terms of the 
running mass. 
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Figure 1. Running of the bottom quark mass 
from low energies to the Mz scale. Upper line is 
the run of rhb{fhb) = 4.39(GeU) with as(Mz) = 
0.112. Bottom line is the run of fh^rhb) = 
4.27{GeV) with a s (M z ) = 0.124. Second pic- 
ture is the difference of both, our estimate for the 
propagated error. 



2. THREE JETS OBSERVABLES AT LO 

Quark masses can be neglected for many ob- 
servables at LEP because usually they appear as 
the ratio m 2 /M^. For the heaviest quark pro- 
duced at LEP, the bottom quark, this means a 
correction of 3 per mil for a quark mass of 5 
(GeV). Even if the coefficient in front is 10 we 
get at most a 3% effect, 1% if we use the bot- 
tom quark running mass at Mz- This argument 
is true for total cross section. However, jet cross 



Table 1 

The jet-clustering algorithms 



Algorithm 


Resolution 


EM 


2 (Pi -Pj)/ S 


JADE 


2(E i E j )(l - cos 0ij)/s 


E 


ipi+Pjf/s 


DURHAM 


2rmn(£f,.E?)(l-costfy)/s 



sections depend on a new variable, y c , the jet- 
resolution parameter that defines the jet multi- 
plicity. This new variable introduces a new scale 
in the analysis, E c = Mzy/yl, that for small val- 
ues of y c could enhance the effect of the quark 
mass as m 2 /E 2 = (m\/M%)/y c - The high pre- 
cision achieved at LEP makes these effects rele- 
vant. In particular, it has been shown j|] that 
the biggest systematic error in the measurement 
of as(Mz) from &6-production at LEP from the 
ratio of three to two jets comes from the uncer- 
tainties in the estimate of the quark mass effects. 

We are going to study the effect of the bottom 
quark mass in the following ratios of three-jet de- 
cay rates and angular distributions 



r>bd — 



rt(j/ c )/r b 




(3) 
(4) 



where we consider massless the d-quark and d 
is the minimum of the angles formed between 
the gluon jet and the quark and antiquark jets. 
Both observables are normalized to the total de- 
cay rates in order to cancel large weak corrections 
dependent on the top quark mass [fi"l| . 

At LO in the strong coupling constant we must 
compute the amplitudes of the Feynman dia- 
grams depicted in figure |^ plus the interchange 
of a virtual gluon between the quark and anti- 
quark that only contributes to the two-jet decay 
rate. In addition to renormalized UV divergences, 
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Figure 2. Feynman diagrams contributing to the 
three-jets decay rate of Z — > bb at order as- 



IR singularities, either collinear or soft, appear 
because of the presence of massless particles 
like gluons. Bloch-Nordsiek and Kinoshita-Lee- 
Nauemberg theorems Q assure IR divergences 
cancel for inclusive cross section. Technically this 
means, if we use DR to regularize the IR diver- 
gences of the loop diagrams we should express the 
phase space for the tree-level diagrams in arbi- 
trary D-dimensions. The IR singularities cancel 
when we integrate over the full phase space. 

Another delicate question is the problem of 
hadronization. Perturbative QCD gives results 
at the level of partons, quarks and gluons, but 
in nature one observes hadrons, not partons, and 
hadronization can shift the QCD predictions. We 
apply to the parton amplitudes the same jet clus- 
tering algorithms applied experimentally to the 
real observed particles, see table [l| Starting from 
a bunch of particles of momenta pi we calculate, 
for instance, = 2(pi -pj)/s, the scalar product 
of all the possible momenta pairs. If the mini- 
mum is smaller than a fixed y c we combine the 
two involved particles in a new pseudoparticle of 
momentum + pj . The procedure starts again 
until all the ytj are bigger than y c . The num- 
ber of pseudoparticles at the end of the proce- 
dure defines the number of jets. The jet clus- 
tering algorithms automatically define IR finite 
quantities. For the moment, we do not enter in 
the question of which is the best jet clustering 
algorithm although the main criteria followed to 
choose one of them should be based in two re- 
quirements: minimization of higher order correc- 
tions and insensitivity to hadronization. If we re- 
strict to the three-jet decay rate the IR problem 
can be overcome and everything can be calculated 
in four dimensions because the jet clustering algo- 



rithms automatically exclude the IR region from 
the three-body phase space. 

For massless particles and at the lowest order 
the EM @, JADE and E algorithms give the 
same answers. Analytical results for the mass- 
less three-jet fraction exist for both JADE-like j|] 
and DURHAM algorithms. A complete anal- 
ysis for the ratios of three-jet decay rates and 
the angular distributions quoted in Eq. ^ and 
H can be found in 0. For practical purposes a 
parametrization of the result in terms of a power 
series in log y c gives a good description . 




Figure 3. Feynman diagrams contributing to the 
three-jets decay rate of Z — > bb at order a|. Self- 
energies in external legs have not been shown. 
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3. THREE JETS OBSERVABLES AT 
NLO 

The effect of the bottom quark mass has been 
studied experimentally by || on the R^ d ratio. 
As we have seen the running of the bottom quark 
mass from low energies to the Mz scale is quite 
strong. The LO QCD prediction for Rl d does 
not allow us to distinguish which mass we should 
use in the theoretical expressions, either the pole 
mass or the running mass at some scale. The 
computation of the NLO is mandatory if we want 
to extract information about the bottom quark 
mass from LEP data. 

At the NLO we have to calculate the interfer- 
ence of the loop diagrams depicted in figure || 
with the lowest order Feynman diagrams of fig- 
ure ^ plus the square of the tree-level diagrams of 
figure ^. The amplitudes in the massless case were 
calculated by The implementation of the jet 
clustering algorithms was performed by [||. 

The main problem that now we can not avoid 
is the appearance of IR singularities. With mas- 
sive quarks we loose all the quark-gluon collinear 
divergences. The amplitudes behave better in the 
IR region. The disadvantage however is the mass 
itself. We have to perform quite more compli- 
cated loop and phase space integrals. Further- 
more, we still conserve the gluon-gluon collinear 
divergences leading to IR double poles. 

The three-jet decay rate can be written as 



T b 3j = C[g 2 v H v (y c ,r b ) + g 2 A H A (y c ,r b )}, 



(5) 



where r b = m 2 /M 2 , C = M z g 2 / '{c 2 w M^){a s / V) 
is a normalization constant that disappear in the 
ratio and gy and g A are the vector and the axial- 
vector neutral current quark couplings. At tree- 
level and for the bottom quark gy = — 1 + 4s^/3 
and gA = 1- Now we can expand the functions 
Hy(A) m a S and factorize the leading dependence 
on the quark mass as follows 



V(A) = A^(y c )+r b B^ A) (Vc,r b ) (6) 
+ ^{A^{y c )+r b B^ (A) {y c ,r b )), 



where we have taken into account that for mass- 
less quarks vector and axial contributions are 



identical^. 

First steep in the calculation is to show the can- 
cellation of the IR divergences in order to build 
matrix elements free of singularities. It is possible 
to do it analytically. However, we knew from the 
beginning IR divergences should disappear |T^| . 
The challenge is in the calculation of the finite 
parts. This calculation is rather long, complex 
and full of difficulties. Strong cancellations occur 
between different groups of diagrams making dif- 
ficult even a numerical approach. We have taken 
as guide line the massless result of although 
the IR structure of the massive case is completely 
different from the massless one. 













TO 5.9 


















= 0.01 




>- 5.8 
p + 5.7 
< 5.6 


























































5.5 












































5.4 
5.3 
























Jad 


B ycut 


= 0.20 














5.2 






















5.1 


















, , , 















0.002 0.003 

r.=m 2 JMl 



0.002 



0.003 

■.ml/Ml 



Figure 4. NLO vector contribution to the three- 
jet decay rate of Z — > bb for bottom quark masses 
from 1 to 5(GeV) and fixed y c in the JADE algo- 
rithm. Big circle is the massless case. 



In figures | and | we present our preliminary 
result for the vectorial contribution to the O(ctg) 
three-jet decay rate of the Z-boson into bottom 
quarks. We have performed the calculation for 
different values of the bottom quark mass from 
1 to 5 (GeF) for fixed y c . We want to show we 
can recover the massless result ||, depicted as a 
big circle, i.e., in the limit of massless quarks we 
reach the (y c ) function. This is our main test 
to have confidence in our calculation. 

In the JADE algorithm we can see that for big 
values of y c the NLO corrections due to the quark 



3 We do not consider the small 0(a|) triangle anomaly 
JL4| . With our choice of the normalization A^(y c ) = 
AJy c )/2 and A^(y c ) = B(y c )/4, where A(y c ) and B(y c ) 
are defined in H. 
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Figure 5. NLO vector contribution to the three- 
jet decay rate of Z — > && for bottom quark masses 
from 1 to 5(Gey) and fixed y c in the E algorithm. 
Big circle is the massless case. 



mass are very small and below the massless result. 
Notice they increase quite a lot for small values of 
y c and give a positive correction that will produce 
a change in the slope of the LO prediction for R^ d . 
In any case we recover the massless limit and a 
linear parametrization in the quark mass squared 
could provide a good description. 

The E algorithm behaves also linearly in the 
quark mass squared although only for big values 
of y c - Corrections in the E algorithm are always 
very strong. The reason is the following, the res- 
olution parameter for the E algorithm explicitly 
incorporates the quark mass, = (pi +pj)/s, 
i.e., for the same value of y c we are closer to 
the two-jet IR region and the difference from the 
other algorithms is precisely the quark mass. This 
phenomenon already manifest at the LO. The be- 
haviour of the E algorithm is completely different 
from the others for massive quarks. It is difficult 
to believe in the E algorithm as a good prescrip- 
tion for physical applications since mass correc- 
tions as so big. However for the same reason, it 
seems to be the best one for testing massive cal- 
culations. 

4. CONCLUSIONS 

We have presented the first results for the NLO 
strong corrections to the three-jet decay rate of 
the Z-boson into massive quarks. In particular, 
extrapolating our result we have shown we can re- 
cover previous calculations with massless quarks. 



Their application to LEP data, together with the 
already known LO, can provide a new way for de- 
termining the bottom quark mass and to show for 
the first time its running. 
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